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Abstract 

A set of MapleV R5 software routines for solving first order ordinary differen- 
tial equations (lODEs) is presented. The package implements the Prelle-Singer 
Method in its original form (see iQ) plus its extension to include elementary func- 
tions (ELFs)0 on the integrating factor (see ||3|, ^, The package also presents a 
theoretical extension to deal with all lODEs presenting liouvillian solutions (LIS) 

H. 
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^E-mails: lduarte@dft.if.uerj.br, lduarte@dft.if.uerj.br, damota@dft.if.uerj.br and jinisk@dft.if.uerj.b: 
^For the definition of elementary function see 0|. 



PROGRAM SUMMARY 



Title of the software package: PSolver. 
Catalogue number: (supplied by Elsevier) 

Software obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland (see applica- 
tion form in this issue) 

Licensing provisions: none 

Operating systems under which the program has been tested: Linux (RedHat 5.2 and Debian 2.0.34), 
Windows 95, Windows 98. 

Programming languages used: Maple V Release 5. 

Memory required to execute with typical data: 32 Megabytes. 

No. of lines in distributed program, including On-Line Help, etc.: 1370. 

Keywords: Elementary fmictions, PrcUe-Singer, first order differential equations, symbolic computing. 
Nature of mathematical problem 

Solving of first order differential equations via the Prelle-Singer Method. 
Methods of solution 

The method of solution is partially the standard Prelle-Singer Method extended to include the case where 
the lODE presents elementary functions. We have also implemented an extension of our own to deal 
with lODEs which solutions are liouvillian. 

Restrictions concerning the complexity of the problem 

If the integrating factor for the lODE in charge has factors of high degree in the dependent and indepen- 
dent variables and in the ELFs appearing on the lODE then the package may expend much computational 
time to catch it. 

Typical running time 

This depends strongly on the lODE. 

Unusual features of the program 

Apart from being a solver of lODEs, implementing the Prelle -Singer approach, our package is a research 
tool allowing the user to follow each and every step of the procedure, for example, which are the eigen- 
polynomials associated with the lODE (see section ||). Besides that, our package proved to be successful 
in solving many classes of lODEs that were missed by some of the best solvers available. Last but not 
least, our package implements a theoretical extension (for details, see section^ and our paper|^) to the 
original Prelle-Singer approach that enhances its scope allowing for it to tackle all lODEs that present 
LIS. 
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LONG WRITE-UP 



1 Introduction 

The fundamental position of differential equations (DEs) in scientific progress has, over the 
last three centuries, led to a vigorous search for methods to solve them. The overwhelming 
majority of these methods are based on classification of the DE into types for which a 
method of solution is known, which has resulted in a gamut of methods that deal with 
specific classes of DEs. This scene changed somewhat at the end of the 19th century 
when Sophus Lie developed a general method to solve (or at least reduce the order of) 
ordinary differential equations (ODEs) given their symmetry transformations [0, ||, 
Lie's method is very powerful and highly general, but first requires that we find the 
symmetries of the differential equation, which may not be easy to do. Search methods 
have been developed |TD|, |TT| to extract the symmetries of a given ODE, however these 



methods are heuristic and cannot guarantee that, if symmetries exist, they will be found. 

On the other hand in 1983 Prelle and Singer (PS) presented a deductive method 
for solving lODEs that presents a solution in terms of ELFs (SELF) if such a solution 
exists The attractiveness of the PS method lies not only in the fact that it is based 
on a totally different theoretical point of view but also because, if the given lODE has 
a SELF, the method guarantees that this solution will be found (though, in principle it 
can admittedly take an infinite amount of time to do so). 

Despite all the good points about the PS-method it presents a shortcoming: it can 
not guarantee to solve lODEs which solutions are not SELF. In order to improve that 
situation, we have developed an extended PS-method. 

So, apart from implementing the PS-method (for the first time, in Maple, with the 
complete approach and the additional bonus of the afore mentioned research tools), we 
present and implement a theoretical extension to the scope of the method. 

The paper is organized as follows: In section ^, we present a short introduction to 
the PS-approach. In the following section, we introduce the main ideas of our extension 
to the procedure. Next, in section ^ we present the package. In sections ^ and |^ we, 
respectively, present some examples and the performance analysis of the package. Finally, 
we end with the conclusion. 



2 The Prelle-Singer Method 

Despite its usefulness in solving lODEs, the Prelle-Singer procedure is not very well known 
outside mathematical circles, and so we present a brief overview of the main ideas of the 
procedure. 

2.1 Rational lODEs 

Consider the class of lODEs which can be written as 

, ^dy_^ M{x,y) 
^ dx N{x,y) ^ ' 
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where M{x, y) and N{x, y) are polynomials with coefficients in the complex field C. 

In [H, Prelle and Singer proved that, if an elementary first integral of (P exists, it is 
possible to find an integrating factor R with i?" G C for some (possible non-integer) n, 
such that 

^ + ^ = 0. (2) 

OX oy 

The ODE can then be solved by quadrature. From (0) we see that 

^BR ^dN ^3R ^dM ^ 
ox ox oy oy 



Thus 



where 



Now let 



D[R] 


fdN 


dM 
dy 


R 


\ dx 


D 


dx ~^ 


dy' 



(4) 
(5) 



R = Ilfr (6) 

i 

where /j are irreducible polynomials and rii are non-zero rationals. From (|^), we have 

D[R] ^ Dju^fn J:^f^'nMm,^^f^ 

R Yiifk'^ rifc fk 

i f i i fi 

From (I), plus the fact that M and N are polynomials, we conclude that D[R]/ R is a. 
polynomial. Therefore, from (|^), we see that fi\D[fi] (i.e., /, is a divisor of D[fi]). 

We now have a criterion for choosing the possible fi (build all the possible divisors of 
D[fi]) and, by using (|) and (0), we have 

If we manage to solve (P) and thereby find n^, we know the integrating factor for the 
lODE and the problem is reduced to a quadrature. Risch's algorithm [0 can then be 
applied to this quadrature to determine whether a solution exists in terms of elementary 
functions. 



2.2 Algebraic and Transcendental lODEs 

Note that the procedure explained above assumed (in eq.|lD that M{x, y) and N{x, y) 
are polynomials in x and y over C[x, y]. 

In order to allow the PS-method to solve differential equations with transcendental 
or algebraic terms, Shtokhamer |TB| regarded the different transcendental and algebraic 
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terms appearing in M{x,y) and N{x,y) as new variables M2, Wn) to be used, in 
addition to the original variables x and y, in the construction of the polynomials fi 
described above. 

Bearing that in mind, it is easy to notice that eq. ^ will become: 

^ = -(^ + y^^]N(xyu, )-(^ + y^^]M(xyu, ) (9) 

A quick inspection of the equation above reveals that its right-hand-side will not 
generally be polynomial, as it would be necessary to apply the PS-method. 
The non-polynomial terms might have two origins: 

^ and ^ may produce functions that are not in {ui, 



• Further than that, ^ and ^ could be rational instead of polynomial. 

We can address the first question by noticing that the Ui, where i = 1, ...n, present 
on the lODE are, by hypothesis, elementary. Therefore, we can construct a complete 
set {ui, ..,Un, ■■Un+m} (which we will call the basis of functions), by including the new 
functions (generated by differentiation as explained above) until we exhaust the process, 
i.e., no new functions appear. 

The second problem can be fixed by multiplying the whole eq.(|^) by all the denomi- 
nators generated by ^ and 

To illustrate these mystical arguments let us use a simple example. Suppose we have 
the following lODE: 

/ /n(x) + sinjx) 

y = r~\ • ^0 

y{x) 

So, we have 

Ui = ln{x), U2 = sin{x). (11) 
Let us first check if the derivatives introduce new functions: 

§fl = i, £p = ^ = 0, = 0. (12) 

OX X ox oy oy 



From eq.(|T2D, one can see that a new function, cos(x), was introduced. So, we have 

Ui = ln{x), U2 = sin{x), M3 = cos{x). (13) 

Now, we would have to, once again, verify if the derivation will introduce new functions. 
In this case, it is easy to see that no more functions will appear, since 

^ = = (14) 

thus completing the set. 

Therefore, we have solved the first problem mentioned above. From eq.(|12D we can 
see that the derivative of ln{x) generates a denominator [x). Then, to ensure that the 
right-hand-side of eq. (||) is polynomial, we would have to multiply the whole equation by 
X (i.e., by the denominator of ^). 

So, we can summarize the approach described in this section by three steps: 
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(15) 



• Then the procedure follows as in the case of the rational lODEs. 

In this section we presented an overview of the PS procedure. Remember that we can 
only guarantee to solve lODEs which solutions are SELF. In the next section, we will 
present a theoretical extension of the PS approach that will permit us to tackle all lODEs 
with LIS. 



3 Extended PS procedure 

In this section, we are going to briefly introduce the ideas behind a theoretical extension 
that will enhance the scope of applicability of the PS-procedure. We are going to use an 
example to make clear what we are talking about. 

Consider the following example, taken from the book by Kamke| [0(eq. k21): 

dy 

— = — ?/sin(a;) + cos(x), (16) 



which solution is 



g— cos{a;) 



y = sin(x) - — m-- i^'^) 

Note that the solution is not a SELF. So the PS-procedure can not guarantee to get 
it. To understand why, let us then consider the integrating factor for this lODE: 



R 



g— cos(x) 



{y - sin(x))^ 



What can we learn from that? Examining R we can see that it is written as a product 
of ELFs. Indeed, we can claim the R is of the form: 

R = l[hT (19) 

i 

where, in this case, hi = y — sin(x) and /12 = e~™^*-^^ Let us now compare this with 
eq.(^). What is the fundamental difference between them? The point is that, in eq.(|TUp, 
not all hi are polynomials over the algebraic extension of eq . (|T^) . In this example, the 
guilty party is /;.2, thus explaining why the PS-procedure would not work. 

Let us introduce our extended PS-procedure that will allow to deal with all those 
cases. Some of what is about to come will involve the Lie theory[^, |^, for lODEs. 

^ A classical reference on the field, generally used as a testing ground. 
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The reader that is not famihar with this should not despair. We are going to explain (in 
English) what every step means. 
Every lODE of the form 

y' = ^= (20) 

^ dx N{x,y) ^ ' 

where M{x, y) and N[x, y) are polynomials of ELFs, presents an infinite number of sym- 
metries. We can guarantee that at least one is of the form: 

[e = 0, r/ = r?(x,i/)] (21) 

In this language, the integrating factor is: 

1 -1 
(22) 

~ M^- Nr] ~ Nr]{x,y) ^ ' 



So, if the integrating factor i?, for the lODE (pO]), is of the form given by eq.(|T9]), then 
?7(x, y) must also be of the form 

ri{x,y) = T[K\ (23) 

k 

since, by assumption, N is polynomial on the algebraic extension of the lODE (^). So, 
for the lODE ([T6|), using eqs. (p^ and (0), we have: 



r/(x, y) = -e"°'(^) {y - sm{x)f (24) 

Applying the Lie theory one can obtain the most general form for an lODE, with a 
given symmetry. For the symmetry presented on eq.(|2l|), this is: 

T- = ^(^' y) r-T^^dy + y)F{x) (25) 
ax J ri[x,yY 

Examining eq.(|TBp, we see that it does not contain e'^°*'^^^ On the other hand, by 
looking at eqs.(p^ and (pSj), it apparently should contain it! How can e™^^^'^ disappear 
from eq.(p5|)?! 

First of all, from the last term of eq. (pSf) , we can see that any multiplicative factor, 
present on ri{x,y), can be eliminated via the arbitrariness of F{x). To be sure that this 
multiplicative factor does not "survive" in the other term (the one with the integral) it 
can not be a general function of x and y. Actually, one can see (from the general structure 
of the term) that this factor has to be an exponential 0. 

Summarizing the ideas presented on this section, we see that the PS-procedure (for 
both the rational and algebraic-transcendental cases) is based on the fact that some 
quantities are divisors of others (see (||) and (|15|)). With the example we have presented, 
we hope to have left it plain to the reader that one might have factors appearing on R 
that are not present on the lODE and, indeed, that those factors have to be exponentials 
where the arguments are functions of x and y. The question that presents itself is: what 
are those functions? By inspecting eqs. (||) and (p!5| ) we can see that they have to be 
taken from the basis of functions of the corresponding lODE (defined, as explained on 
section |2.2|, by the functions present on and M of the lODE). 
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We would like to point out that the example dealt with on this section (and the 
reasoning following it) is just a part of the whole body of ideas and arguments we present 
on to conclude that there is a possible extension to the PS-procedure. Nevertheless, 
it is a complete part in the sense that it should be enough to "convince" the reader of 
the importance of this extension. In other words, we hope to have showed that, without 
this, a whole class of lODEs would elude the PS-procedure. We will, later on this paper, 
show some more examples only solved by this PS-extended-procedure and comment on 
our computational implementation of these ideas. 



4 The PSolver package 
Summary 

A brief review of the commands of the package is as follows:^ 

• PSolve solves lODEs, using the Prelle-Singer approach. It deals not only with ra- 
tional lODEs but also with the ones containing ELFs (algebraic and transcendental 
functions). 

• PSolvexp extends the applicability of the PS-procedure to all lODEs with LIS. 

• Intf act returns an integrating factor for the lODE. 

• Dop constructs the V operator associated to the lODE. 

• Basis constructs a list of equations of transformation which name a set of algebraic 
and transcendental functions (e.g., [ui = fi{x,y), ...]). Those functions are, see 



sections |2.2| and [4.3| , what we call the basis of functions for the corresponding lODE. 



• dBasis presents the list with the partial derivatives of the basis of functions (e.g., 
[dui/dx = -U2, ...]). 

• Eigenpol calculates all the polynomials fi that are divisors of T>fi, i.e., that satisfy 
the equation Vfi = gifi. 

• EigenPval calculates the polynomials gi that are the "eigenvalues" for the equation 
T^fi = 9ifi- 

Description 

A complete description of the PSolver package's commands is found in the on-line help. 
Here we present the most relevant part of that description. 



■^This subsection and the next one may contain some information already presented in the previous 
sections; this is necessary to produce a self-contained description of the package. 
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4.1 Command name: PSolve 



Feature: This command applies the PS-procedure to solve lODEs. 
Calling sequenc^: 

> PSolve(lode,n,optional_parameters) ; 

Parameters: 

lode - a first order ordinary differential equation. 

n - an integer which sets the highest degree the command will consider for finding the 
polynomials /j. 

Optional Parameters: 

NumberF=integer - Tells the program the number of fi it is going to 
use to look for the integrating factor. 



NumberF=All - Tells the program to use all fi in searching 

for the integrating factor. 



AllSolutions 



- The absence of this argument truncates the searching for all 
possible fi in a way that, for the majority of cases, still allows 
the command to find the integrating factor but with considerable 
gain in time. 



Extended 



- Indicates to the program that it should use the set of routines 
that look for the exponentials in order to calculate the integrating 
factor. 



Guess - It activates an algorithm for fast determination of fi 

candidates. 

Synopsis: 

The PSolve command is a part of the PSolver package, which, to brutally summarize 
things, is designed to solve lODEs with LIS. The command implements the PS-procedure 
and an extension of our own (see section ^ which can be divided into five steps: 

• First, from the lODE, one can determine the basis of functions. 

• In second place, from the basis (mentioned above) one constructs the V operator, 
see sectionp^. 

• Following that, we can calculate the "eigen-polynomials" and "eigenvalues" (which 
are themselves polynomials) for the operator V. 

• From these results, we find (or not, if it does not exists) the integrating factor for 
the lODE. 

^In what follows, the input can be recognized by the Maple prompt >. 
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• Finally, one can then use the integrating factor to find the desired solution for the 
lODE. 

The arguments 

The first argument of PSolve is the lODE. The second argument is an integer that 
tells the command the highest degree it is going to use in the search for the integrating 
factor. We would like to remind the reader that the PS-procedure is a semi-decision 
one|Q, 1^, in other words, if the integrating factor has not be found for this degree, it does 
not exists at this level. But, it could as well exist for higher degree, we have to keep 
looking. 

The optional arguments 

The optional arguments can be given alone or in conjunction and in any order. In 
order to improve its practicability (speed combined with efficiency) PSolve, by default, 
uses seven to construct the integrating factor. Why is that so? Well, this process of 
calculate the integrating factor can be broadly divided into two steps: first one has to find 
all possible /«, then use them to look for the integrating factor. This process of looking for 
the IF is "computationally expensive" and, furthermore, we have noticed that it is seldom 
the case where all the fi are needed to do the job. We have also noticed that, generally, 
using the "simpler" fi was enough. What do we mean by "simpler"? We call a simple 
fi the one which leads to a faster calculation. Luckily for us, there is a direct relation 
between this desired property and the form of the fi. Actually, we have observed that 
the "smaller" the fi the faster is the computation time. So we introduced a algorithm for 
ordering the possible fi by "size"Q The default number seven was defined heuristically. 
Of course there will be the odd case where seven won't be enough, so we let the door 
open for the user to increase this number via the optional argument NumberF=integer. 

How about the other step? Can we improve the efficiency of our package in the finding 
of all /j? In fact, for the lODEs which involve radicals, the problem of computational- 
time-expenditure can be even greater. For instance, in a particular example, the program 
took several minutes to find some forty fi. It turned out that we developed a different 
algorithm that, although loosing some possible solutions, still produces enough of them 
(for a great number of cases) thus allowing the program to find the IF. For the remaining 
(few) cases, the user may then apply the optional argument AllSolutions which will tell 
the program to use the first (slower) algorithm. 

The importance of the option Extended lies in the fact it implements a body of new 
ideas making it possible for the user to extend the PS-procedure to a greater range of 
lODEs, to cover all the ones with LIS (not only with SELFs). See section ^. 

Finally, the Guess option implements an algorithm for fast pattern recognition de- 
termining possible fi candidates of high degree, thus shortening the computational time 
required by the PS-procedure. This algorithm will help in some (difficult) cases but it 
complicates the general one and so it is optional. 

4.2 Command name: Intfact 

Feature: Calculates a integrating factor (R) for the lODE. 

^This definition of "size" and the creation of the routines to measure it are a part of our package. 
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Calling sequence: 

> Intf act (lode ,n, optional_parameter) ; 

Parameters: 

The parameters lode and n have the same meaning as explained aboveQ. 

Optional Parameter: 

NumberF=integer - Tells the program the number of fi it is going to 
use to look for the integrating factor. 



NumberF=All - Tells the program to use all fi in searching 

for the integrating factor. 



AllSolutions 



- The absence of this argument truncates the searching for all 
possible fi in a way that, for the majority of cases, still allows 
the command to find the integrating factor but with considerable 
gain in time. 



Extended 



- Indicates to the program that it should use the set of routines 
that look for the exponentials in order to calculate the integrating 
factor. 



Guess - It activates an algorithm for fast determination of fi 

candidates. 

Synopsis: 

Actually this command embodies the heart of the package and that is so for, basically, 
two reasons: The brilliance of the PS-procedure lies exactly on the calculation of the 
integrating factor. After that step, the finding of the solution is no novelty. The second 
reason is that our extension to the PS-procedure can be also found mostly on this com- 
mand (via the optional parameter mentioned above). Again, after we find the integrating 
factor, it is straightforward to look for the solution to the lODE. 



4.3 Command name: Basis 

Feature: Determines the basis of functions for the lODE. 
Calling sequence: 
> Basis(lode) ; 
Synopsis: 

With the somewhat more specialized user in mind, we implement a set of commands to 
help the research with lODEs. The first one is this command that generates, from the 
ELFs present on the lODE, a set of functions that is complete under differentiation. 

^In what follows, these parameters will always have this meaning. 
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4.4 Command name: DBasis 

Feature: Returns a list with the partial derivatives of the functions in the Basis. 
Calling sequence: 

> DBasis (lode) ; 

Synopsis: 

This is another one of the commands to help in the research on the structure of lODEs. 
It generates a table with all the first partial derivatives of the Basis, written in terms of 
this Basis and of the independent and dependent variables. This will be used to construct 
the V operator. 

4.5 Command name: Dop 

Feature: Returns the V operator (for the lODE) that is a essential ingredient in applying 
the PS-procedure. 

Calling sequence: 

> Dop (lode) ; 

Synopsis: 

This command returns the T> operator associated with the lODE. The operator is 
written in the basis of functions for the lODE (see section |]). 

4.6 Command name: Eigenpol 

Feature: Returns the "eigen-polynomials" for the T> operator. 
Calling sequence: 

> EigenpoKlode ,n) ; 

Synopsis: 

This command finds all the polynomials fi (of maximum degree n) satisfying the 
equation Vfi = gifi and those will be the candidates to build the integrating factor for 
the lODE. 

4.7 Command name: EigenPval 

Feature: Returns the "eigenvalues" for the V operator. 
Calling sequence: 

> EigenPvaldode ,n) ; 
Synopsis: 

This command finds the polynomials gi that play the role of eigenvalues for the equa- 
tion Vfi = gifi. 
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5 Examples 



In this section, wc arc going to exemplify some aspects of the routines in the package 
and show that, besides being vahd only on the basis of being an implementation of the 
regular PS-procedure (and that this fact already greatly enhances the solving of lODEs 
capabilities within the Maple environment), our program proved to be a successful attempt 
of implementing news ideas, extending even further the applicability of the PS-procedure. 

In order to avoid boring the reader to death and to allow him/her to focus on the 
relevant aspects, we are going to concentrate (try to) in a few lODEs that will serve as 
examples for as many features as possible. 

5.1 Usage of the package commands 

Let us use a simple lODE to show the commands in action. Consider the following lODE: 

dy ^ y (cos(x) + ye"^ + 1) 
dx cos{x) 

Usually, people only want to find the solution for the lODE. So, that would simply 
require the input lines: 

> eq := diff(y(x),x) = y(x)*(cos(x)+y(x)*exp(-x)+l)/cos(x) : 

> PSolve(eq,l) ; 

]n{y{x) + e") - ln(tan(l/2 x) - 1) + ln(tan(l/2x) + 1) - \n{y{x)) - _C1 = (27) 

Sometimes the user is interested on having a look at just the integrating factor, so 
he/she could type: 

> R := Intf act (eq, 1) ; 

For the 1 ODE in the form 

d , X y(a:) (cos(x)e^ + yix) + e^) 

—y\x) = . . — — 

dx cos(x)e^ 

the integrating factor will be 



y cos(x) {y + e^) 

Please notice that the program informs the user the format in which it is considering the 
lODE to be in, that is important since it may be the case where the input lODE gets 
transformed by the simplification routines and is regarded in a different shape (of course 
it is still the same lODE) and that will affect the integrating factor. 

The package can also provide information about the basis of functions (that can be 
useful in analyzing the structure of the lODE). 

> basis := Basis (eq); 
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basis = [ui = cos(x),M2 = sin(x),'U3 = e^] 

For the user that is interested in the hmits of apphcabihty of the PS- procedure or that is 
otherwise involved in more pure research in the algebraic properties of lODEs, that can 
prove to be a very valuable piece of information. The same applies to the commands that 
follow. 

For instance, the following command to be exemplified [dBasis) provides essential 
information in order to apply the PS-procedure (please refer to the above definition of 
basis and, for further detail and meaning, see section p^ : 

> d_basis := dBasis (eq); 

d d d 

dMsis = [—ui{x, y) = -U2, ^"2(3;, y) = ui, ^^3(3;, y) = M3, 

Odd 
— y) = 0, ^"2(3;, y) = 0, ^^3(3;, y) = 0] 

The above result is used by the following command which defines the V operator. 

> D_operator:= Dop(eq); 

^ d . , (9 d 

D_operator = w uiu^—w + y {U1U3 + y + U3) t— w — U2U1U3— — w+ 

ox oy oui 

Ul U3-—W + U3 Ui^—W 

OU2 OU3 

The user who decides to use this command has to bear in mind that it returns a procedure. 
Ops, what does that mean? It defines a mapping to be used later on by the user. To 
clarify what we mean, the best to do is to see it in action in a simple case. Let us apply 
the recently defined D_operator to (for example) Ui. One would get: 

> example := D_operator(u[l] ) ; 

example = —U2U1U3 

Of course, one has to apply Djoperator to quantities written on the basis, i.e., as a 
function (for our example) of x, y, Ui, U2 and M3. 

The next two commands are very specialized (for the user really into lODE research). 
Please see sections ^ and ^ for recalling the meaning of them. For our particular example, 
we have: 

> eigen_p := Eigenpol (eq, 1) ; 

eigenjp = [y, e^, cos(x), y + e^, cos(x) — /sin(a;)] 

> eigen_p_val := EigenPval(eq, 1) ; 
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eigenjpjval = [cos(x)e^' + y + e^, cos(x)e^', —e^ sm{x),y + cos(x)e^', — / cos(x)e^] 

Let us have a closer look at this results. The two lists above are consistent in the sense 
that their order is related, i.e., the "eigen-value" associated with y is cos(a;)e^ + y + e^, 
the "eigen-value" associated with is cos(a;)e^ and so on. In the language of the basis 
of functions, that would translate to, for example: 

T^[y] =yiuiU3 + y + u3) 

As previously mentioned, we have used a single lODE to exemplify the usage of 
the package commands, we hope to have clarify that usage and helped the reader to 
understand further the PS-procedure. Now we are going to further exemplify the package 
but with a different approach. We will show that our package complements the solving 
capabilities of the existing solvers in Maple. We will divide that into two main groups: 
the first talks about the fact that our implementation solves lODEs, with SELfs, that are 
missed by the powerful solvers on the MAPLE algebraic package. The second one, deals 
with our theoretical extension to the PS-procedure (and our implementation of that) and 
tackles lODEs with LIS. 



5.2 lODEs with SELFs 

The solver (in MAPLE) is the dsolve command and it is consider by many the most 
powerful solver of ordinary differential equations commercially available. It fails to solve 
the following two examples: 

egi = ^ = ^ + 9^-6ye^ + y^ (28) 
ax 

dy (y^ ln(x)^ + 4?/ln(x)^ + 4 ln(a;) + y^) 
""^^ ~ rf^ ~ {y\n{xY + 2fx 

Our program finds the following integrating factors for the above lODEs: 

1 

R2 = — 7 

xy^ 

which, in turn, lead to the following solutions: 

soh = -l + 9xe'' -Sxy^ + dCie"" -3 _C1 y^ = 

soh = -y^ \n{xY -6y^ H^Y -I2y ln{xf - 6 \n{x)y^ -8 + 6 Cgy^ = 

It is worth reminding the reader that the dsolve command brings within it many 
algorithms implementing the Lie Method for solving ordinary differential equations 



TT|. But, for an infinity of cases (of which the two equations above are an example), the 
symmetries can not be found by the this heuristic approach. On the other hand, the 
PS-procedure, being of a semi-decision nature, guarantees that, if the solutions exists it 
will be caught (it is only a matter of computational power). 
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5.3 lODEs with LIS 



In this subsection, we are going to exemplify the usage of the theoretical extension to the 
PS-procedure we have explained on section which allows the procedure to be extended 
to deal with lODEs with LIS. 
Here are two such examples: 

dy y + l + e^x^ 

eq2 = — = ^ (30) 

dx x'^y 

dy — cos{x)xy — cos(x)a; — cos(x) + y + 1 + xe^ 
dx 1 + xy 

Our program finds the following integrating factors for the above lODEs: 

1 



Ra 



X e^e^ 
1 



which, in turn, lead to the following solutions: 

soh = -e^+""' C3 + l + y + e^+""' J e"+""'rfa; 

sok = {-xy -1-x) e-y-''''^''^ - J xe-^'^^^^^rfx - Q = 

As previously mentioned, these equations would elude the regular PS-procedure. Fur- 
thermore, these examples (again they exemplify an infinitude of such cases) also escape 
solution via the dsolve command. So, the implementation of the extended PS-procedure 
0, besides enlarging even more the scope of soluble lODEs, now allows for the fact that 
all lODEs with LIS (which include the ones with SELFs) can be treated in the context 
of a semi-decision procedure. 



6 Performance 

In this section, in order to show the efficiency of the ideas hitherto exposed, we are going 
to use the equations on the book by Kamke [Q. Which equations should we choose 



from the book? There are 367 lODEs of first degree on the book. In principle, this 
would be our testing arena. However, among those equations there are, for example, 
some quadrature, some exact lODEs, etc. So, we see no point in using those in our test, 
therefore we will exclude them (see table |I]). There is another important restriction one 
should consider before embarking on the test, some lODEs are "out-of-scope" , i.e., they 
can not be solved by the PS-procedure. More specifically, we mean lODEs with arbitrary 
functions, with solutions in terms of special functions and some lODEs with LIS (see 
table m). The latter case will be only excluded from our test for the usual PS-procedure. 
Surely, we will include those in our test for the Extended-PS (See section |6^ ). 
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Quadrature, Exact 


1, 4, 7, 12, 17, 23, 26, 39, 57, 59, 61, 69, 70, 71, 73, 76, 89, 
90 131 150 154 200 209 223 227 229 245 248 251 
263, 267, 270, 271, 273, 274, 285, 288, 289, 290, 298, 299, 
300, 305, 307, 308, 309, 310, 311, 322, 330, 335, 336, 341, 
348, 352, 353, 355, 356, 358, 360, 361, 366 


ODEs presenting 
Arbitrary Functions 


10, 11, 16, 33, 34, 35, 49, 50, 51, 53, 54, 55, 56, 72, 74, 79, 80, 
84, 85, 86, 87, 110, 126, 127, 128, 201, 202, 212, 219, 230, 250, 
268, 269, 330, 331, 365, 367 


ODEs with Special 
Functions in the 
Solutions 


13, 14, 24, 25, 30, 36, 37, 40, 43, 45, 47, 48, 63, 66, 82, 88, 95, 
99, 100, 105, 107, 111, 121, 139, 144, 145, 146, 147, 157, 166, 
168, 169, 176, 179, 203, 205, 206, 234, 237, 253, 265 


ODEs with LIS 
out of the scope 
of the PS-method^ 


5,18, 20, 21, 22, 27, 28, 46, 83, 129, 133, 164, 235, 343, 351 



Table 1: Classification of ODEs (in Kamke) that will not be considered on the procedure 
testing (the last row will be considered on the testing of the Extended PS-procedure). 



6.1 Testing the usual PS-procedure 

From table |I| we see that our arena for testing the usual PS-procedure consists of 215 
lODEs. Below, see table |^, we present the results of the application of our package to 
those lODEs. 

To complete our test we have to mention the lODEs not present on table Those are 
solved using the optional arguments described on section^. The lODEs 62, 112, 113, 114, 
115, 116, 333, 357 are caught using the optional parameter AllSolutions (115, 333 need, 
besides that, option NumberF=All). The lODEs 151, 182, 185, 211, 257, 327 are caught 
by the single use of NumberF=All. The lODEs 38, 42, 44, 77, 78, 106 are solved using 
the special algorithm that allows for fast solution of "hard" cases (described on section 

This algorithm is enabled by the use of the optional parameter. Guess. 

Finally, the lODEs 52, 81, 142, 173, 184, 186-189, 266, 292 exceed the allowed time 
in our testing (500 seconds). 

6.2 Testing the Extended PS-procedure 

From table 0, we come to the conclusion that 15 lODEs should be considered for our test- 
ing procedure. These are lODEs that would escape solution via the usual PS-procedure. 
The lODEs are numbered 5, 18, 20, 21, 22, 27, 28, 46, 83, 129, 133, 164, 235, 343, 351. 
From those, all but two are solved by our implementation. The lODEs 22, 46 exceed the 
allowed testing time of 500 seconds. 

To conclude our section on the performance of the package, let us just comment that, 
apart from the examples we gave on section ^, which are representatives of an infinitude 
of cases that are missed by the powerful solver in Maple (dsolve), the lODEs numbered 
42, 151, 152, 185, 257, 350, 351 are also solved by our package and missed by dsolve. 
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Average Time 
of solution 


ODEs 


N = 1 


1 3 s 


2, 3, 6, 8, 9, 19, 29, 31, 58, 60, 64, 65, 67, 68, 75, 
91, 92, 93, 94, 96, 97, 98, 101, 102, 103, 117, 118, 119, 120, 
122, 123, 124, 125, 130, 132, 134, 135, 136, 137, 138, 148, 
149, 152, 153, 155, 156, 158, 159, 160, 161, 162, 165, 167, 
171, 174, 175, 177, 178, 180, 183, 190, 191, 192, 193, 194, 
197, 198, 199, 204, 207, 210, 213, 214, 215, 216, 218, 221, 
222 224 225 226 228 231 232 233 236 238 239 240 

241, 242, 243, 244, 246, 247, 249, 252,254, 255, 256, 
258 259 260 261 262 264 272 275 276 277 279 280 
281, 282, 283, 284, 286, 287, 291, 293, 294, 295, 296, 297, 

301, 302, 303, 306, 313, 314, 315, 317, 318, 319, 321, 
323, 324, 325, 326, 327, 328, 329, 332, 333, 334, 338, 339, 
342, 344, 345, 346, 347, 349, 350, 354, 357, 359, 362, 363, 364 


N = 2 


3.3 s 


15, 41, 104, 108, 109, 140, 141, 143, 163, 170, 195, 208, 
217, 220, 304, 312, 340 


N = 3 


26.5 s 


172, 181, 278, 320 



Table 2: This table presents the average time it took our program to solve lODEs in the 
book by Kamke. The table is organized by the value of N, the degree of the /j used to 
construct the integrating factor. 



Besides that, the implementation of the theoretical extension allows us to substantially 
decrease (in some cases) the time needed to solve the lODE. For instance, from the list 
of lODEs that exceed the 500 seconds time limit we have set, the lODE 142 is solved by 
the extended procedure in only 3.4 seconds. 

7 Conclusions 

In this paper, we presented an implementation of the PS-procedure. As mentioned before, 
this method has the advantage of being of a semi-decision nature, i.e., it guarantees that, 
if the solutions exists it will be caught (it is only a question of time). Furthermore, we have 
made available to the user all the intermediary steps of the PS-procedure: the integrating 
factor, the V operator, the basis of functions for the lODE etc (see section H). 

We have also implemented a theoretical extension we have developed^. This enlarges 
the scope of PS-procedure to allow it to tackle all lODEs with LIS (that contain the set of 
the lODEs with SELfs previously tackled by the usual PS-procedure). As an additional 
bonus to all this, we have enhanced the lODE solving capabilities in Maple, since this 
implementation of our extension enables the package to solve lODEs previously missed 
by the powerful solver in Maple (dsolve). 

As further development for this work, we are currently working on extensions to the 
package to tackle lODEs containing arbitrary functions and on the implementation of the 
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Collins approach (and references therein). 
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